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Tight constraints on F- and D-term hybrid inflation scenarios 
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We use present cosmological data from the cosmic microwave background, large-scale structure 
and deuterium at high redshifts to constrain supersymmetric F- and D-term hybrid inflation sce- 
narios including possible contributions to the CMB anisotropics from cosmic strings. Using two 
different realizations of the cosmic string spectrum, we find that the minimal version of the D-term 
model is ruled out at high significance. F-term models are also in tension with the data. We also 
discuss possible non-minimal variants of the models. 
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Introduction: Present observations of the Universe pro- 
vide strong evidence for a topology which is close to spa- 
tially fiat and a close to scale invariant spectrum of den- 
sity fluctuations. These are both predications of the sim- 
plest models of inflation based on the slow-roll approx- 
imation and hence provide prima-facie evidence for an 
inflationary epoch in the Early Universe. One of the ba- 
sic features of inflation is that the potential needs to be 
very shallow in order to give rise to perturbations with 
the observed amplitude. This can be difficult to arrange 
in many particle physics motivated scenarios since it typ- 
ically requires parametric fine-tuning in the theory or a 
trans-Planckian field-range for the inflaton. 

One class of models which naturally circumvents this 
issue are the Supersymmetric (SUSY) Hybrid Inflation 
scenarios which have F- and D-term variants In 
these models the tree level potential is completely flat 
and the slow-roll is induced by the radiative corrections. 
There are two key observational features of the minimal 
version of these scenarios: a spectral index of density 
fluctuations, ns, which is 0.98-1.0, and cosmic strings are 
often produced during the phase transition which ends 
inflation for example, if the broken symmetry is U{1). 
For typical model parameters the amplitude of cosmic 
microwave background (CMB) anisotropics from strings 
will be ~ 10% of those produced by inflation 

This is particularly interesting since the data from the 
Wilkinson Microwave Anisotropy Probe (WMAP) 3 year 
release suggested that ns = 0.958 ± 0.016 which ap- 
pears to disfavor such models Q . In subsequent work we 
pointed out that this discrepancy could be ameliorated 
by the inclusion of a component of the CMB anisotropics 
from cosmic strings 0]. In particular, we found that the 
minimal version of the F- and D-term scenarios, which 
have 6 free parameters, fitted the data as well as or better 
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than the standard 6 parameter fit. Similar results, based 
on a different cosmic string spectrum, were reported by 
Bevis et al, although they did not specifically refer to the 
SUSY hybrid inflation models, rather they investigated a 
6 parameter model with the string tension, G/i, allowed 
to vary and ns = 1 fixed Q. 

Since then cosmological data has moved on. WMAP 
have released their 5 year data @, other CMB data 
on smaller angular scales has been published [l(|, the 
SDSS has been used to compute the matter power spec- 
trum [ll| and tight constraints have been derived on 
the baryon density, quantified by the value relative to 
the critical density, ilb, using Big Bang Nucleosynthe- 
sis (BBN) from measurements of deuterium at high red- 
shift [HI ■ The last two of these are particularly interest- 
ing in this context since we found in our earlier analysis 
that models which fitted the data due to the inclusion a 
string component with ns ~ 1 typically have values of 
h = Ho/100 kmsec -1 Mpc -1 and Sl^h 2 which arc larger 
than the standard values 

In this letter we will revisit the constraints on the 
SUSY hybrid inflation models focusing on the minimal 
F- and D-term variants and on those that are extended 
by a possible additional mass-square term for the infla- 
ton |0 - [l6l |. We will take into account the cosmic string 
component using two possible models. In model A we 
use the unconnected segment model [l7| (USM) to pro- 
duce the string spectrum using an evolution model for the 
correlation length, rms string velocity and string wiggli- 
ness. This spectrum was used in our previous work [7| 
and also in refs. [l8|. We also include results based on 
another spectrum created using the USM, which we will 
call model B. This spectrum, with different parameters 
from model A, was designed to match the results of the 
Abelian-Higgs simulations reported in ref. Q (see ref. [lj| 
for a discussion) . It is our opinion that model A is likely 
to be closest to the true cosmic string spectrum, but the 
two spectra span the different possibilities reported in the 
literature. We note that since the cosmic string spectrum 
will be ~ 10% of the total, the level of accuracy required 
is much less than for the adiabatic spectrum, and the two 
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models give rise to qualitatively similar results. 

Methodology: Our basic methods will be essentially the 
same as those described in rcf. Q. The F-tcrm super- 
potential is given by W = kS(GG — M 2 ) where S is 

the singlet inflaton superfield and G, G are the waterfall 
fields responsible for the topological defects Q. We in- 
clude the radiative corrections and minimal supergravity 
corrections, but ignore the curvature and tadpole terms. 
These terms were shown to have negligible effect in pre- 
vious work, where we assumed a positive sign for the tad- 
pole term. For the D-term model we use W = kSGG. 
Spontaneous breaking of U(l) is induced by the D-term 
contribution g 2 {\G\ 2 - \G\ 2 + M 2 ) 2 /8 to the potential. 
Without supergravity corrections, the predictions for in- 
flation and strings from the D-term model do not depend 
on the gauge coupling g. Taking account of minimal su- 
pergravity corrections, we find that for values of g > 0.1, 
the spectrum receives an additional blue-tilt, which cor- 
responds to a less-good fit to the data Q. When com- 
pared to small g, the hybrid inflation prediction for ns 
lies substantially above the central value from the stan- 
dard 6-parameter fit. In this work, we neglect minimal 
supergravity corrections, corresponding to a restriction 
of parameter space to g ^ 0.1. Therefore in both the 
F- and D-term cases, the inflation model is specified by 
two parameters: the dimensionless coupling constant, k, 
and the mass scale, M. It is possible to compute the 
three observable quantities rig, G/i and the scalar ampli- 
tude Vb, from these two parameters, and hence the fit 
is for 6 parameters (with the baryon density, fl^h 2 , cold 
dark matter density fl c h 2 , acoustic scale, 9a and optical 
depth, t, being the other 4). The reheat temperature is 
fixed to be Tr = 10 9 GeV but allowing this to vary only 
has a minor effect on our results. 

We refer to more details on the models to ref . Q , but in 
the present context, it is useful to recall the basic features 
distinguishing the F- and D-term models phenomcnolog- 
ically. For n 2 / g 2 < 1/2 the string tension for the same 
value of M is larger in the D-term model than the F- 
term. Hence, in order to comply with upper bounds on 
the string contribution to the CMB, D-term models pre- 
dict a smaller M, which also leads to a close to scale 
invariant spectrum, ns w 1. 

We use COSMOMC [H to perform the analysis. 
We include data from WMAP5 other CMB exper- 
iments [l(| and the SDSS [ll|. In addition we will also 
include constraints on Q\,h 2 from metal poor damped 
Lyman-a systems. These have been used to compute 
the deuterium abundance at z sa 2.6. One finds that 
\og(D/H) = -4.55 ± 0.03, which implies that n h h 2 = 
0.0213 ±0.0010 [3. 

Results: We find that the 4 cosmological parameters 
are close to the values preferred by the standard 6 param- 
eter fit. In Fig. [TJ we present 2D likelihoods for various 
parameters of the F- and D-term models for cosmic string 
spectra A and B. The best-fitting values of M and n are 
given in table [TJ We see that the values are not the same 
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FIG. 1: 2D likelihoods for D-Term (red, smaller values of k) 
and F-term (blue, larger values of k) for string model A (left 
set of two) and B (right set of two) using CMB+SDSS data. 



for the two string spectra. The central values for model 
A have changed from those reported in ref. Q due to the 
inclusion of the new CMB and the SDSS data, although 
they remain within ~ la. The size of the errorbars have 
not changed significantly. The results for model B are 
somewhat different reflecting different ranges of G/z and 
ns which are favoured using the different string spec- 
tra [lj|. Of particular importance is that model A re- 
stricts ns = 0.958 ± 0.013 in a 7 parameter fit including 
G/j, as the extra parameter, which is very similar to the 6 
parameter case. However, ns = 0.970 ± 0.016 for model 
B, which is more compatible with ns = 1. 

One can define the relative likelihood of two models, 
labelled 1 and 2, by A x 2 = -2 log(£i/£ 2 ). We will take 
model 1 to be the standard 6 parameter model using 
ns and Aq to define the initial power spectrum. Pos- 
itive values of this quantity indicate that model 2 fits 
the data better than the standard, whereas negative val- 
ues indicate a less good fit. This quantity was used in 
ref. @ to quantify the level to which various variants 
of the standard cosmological model (for example, those 
with no reionization, or no dark matter) were poor fits to 
the data, We present the computed values of A% 2 for the 
F- and D-term models using the two cosmic string spec- 
tra in table |TTJ We see that if only CMB data is included 
there is little to tell between the 6 parameter model and 
the F/D-term models, although the D-term model using 
string spectrum A is the least good fit. If one includes 
more data from SDSS and BBN, typically the fit of the 
hybrid inflation models become less good. 

Concentrating on the case of CMB+SDSS+BBN, it is 
clear that the D-term model using both string spectra 
is a much worse fit than the 6 parameter model with 
A\ 2 — —17.6 for string spectrum A and -12.2 for B. This 
shows that the D-term model is severely in tension with 
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logio « 


A//10 16 GeV 


F(A) 


-2.79 ±0.30 


0.450 ± 0.062 


F(B) 


-1.94 ±0.34 


0.563 ± 0.038 


D (A) 


-4.55 ±0.21 


0.196 ±0.030 


D(B) 


-3.84 ±0.23 


0.330 ± 0.047 



TABLE I: Constraints on minimal F and D-term parameters 
using string models A and B. 

the data. This is because the best-fitting D-term models 
and those with ns = 1, in the absence of the BBN prior, 
typically favour values of fl h h 2 ~ 0.024 - 0.025. The 
situation is more complicated in the F-term case. In both 
cases the F-term model is disfavoured relative to the 6 
parameter model, but for string spectrum B the A\ 2 is 
only —5.6. This would not be enough to exclude the 
model with high confidence, whereas for spectrum A the 
A\ 2 — —12.2, which would. We have already alluded to 
the reason for this: spectrum B allows for larger values 
of rig than the spectrum A. 

In order to calibrate our interpretation of these values 
we have done two things. First, we computed Ax 2 for 
a 5 parameter model with r = 0. For each combination 
of data, we find Ax 2 « —20, with roughly 5er signifi- 
cance that r / (for example, CMB data alone gives 
r = 0.090 ± 0.017). We have also done the same for 
models with ns = 1, that is a Harrison-Zcl'dovich (HZ) 
spectrum. For CMB±SDSS data we find A X 2 = -11.8, 
corresponding to the HZ model being excluded at 3.4cr, 
and for CMB±SDSS±BBN data A X 2 = -16.8, corre- 
sponding to 4.1cr (ns = 0.951 ± 0.012). These values are 
presented in table [TTJ 

Comparing the A\ 2 for the F- and D-term models we 
can conclude that the D-term model is excluded at a 
level somewhere in the range of 3-4er even if strings are 
included, with the precise level dependent on which string 
spectrum is correct. Meanwhile the F-term model would 
be excluded at around 3<r if string spectrum A is correct, 
but cannot be excluded is spectrum B is correct. 

The statistical strength of our results is very strongly 
related to the value of D/H measured in ref. [l2| and 
the subsequent interpretation in terms of BBN. In order 
to explore the effects of systematic offsets in the value of 
fl^h 2 on our results we have computed the value of A% 2 
for a prior of cjb = Qbh 2 = ± 0.001 (i.e. a different 
central value with the same statistical error as ref. [l2l]). 
For the D-term model with string spectrum A, the vari- 
ation in x 2 for CMB±SDSS±BBN can be approximated 
by Ax 2 ~ 3160 Wb — 84.5. Even if the central value were 
to increase by 0.002, the model would still be in tension 
with the data, Ax 2 w — 11. 

Non-minimal models: Non-minimal F- and D-term 
models can be considered; these typically involve the in- 
clusion of a negative mass-square term in the potential 
and can allow for values of ns lower than 0.98 due to an 
enhanced negative curvature of the potential. In addi- 
tion, at the point of horizon exit, the potential can be 
effectively flattened. This leads to a lower energy-scale 





CMB 


+SDSS 


+BBN 


±SDSS/BBN 


r = 


-21.7 


-20.4 


-20.8 


-20.2 


HZ 


-9.4 


-11.8 


-14.6 


-16.8 


F(A) 


-1.2 


-4.4 


-10.2 


-13.2 


F(B) 


1.2 


0.2 


-4.0 


-5.6 


D(A) 


-4.0 


-8.6 


-14.2 


-17.6 


D(B) 


0.2 


-2.8 


-9.6 


-12.2 



TABLE II: Ax 2 = -21og(£i/£ 2 ) for minimal 6 parameter F 
and D-term scenarios. These are defined as model 2; model 
1 is the standard 6 parameter power law. We also present 
results relative to the 5 parameter r = and HZ models. 

of inflation and therefore to lower values for M and a 
smaller string tension. In ref. Ml we considered a par- 
ticular class of F-term models [13| characterized by the 
addition of the potential term c 2 H H 2 S 2 , where H denotes 
the Hubble rate. This "Hubblc-dcpendent" mass term 
generically arises in the presence of a non-minimal Kahler 
potential from the F-term that drives inflation. 

In the absence of F-terms, even for a non-minimal 
Kahler potential, no such additional mass-square correc- 
tion can occur. This is why for the simplest realizations 
of D-term inflation, the exclusion of mass-square correc- 
tions is well-motivated. However, there are particular 
models where sub-dominant F-terms {e.g. from compact- 
ification in string models or from sneutrino potentials), 
that do not drive inflation, may be present [l4T - fl6| . These 
can lead to mass-square corrections for the D-term po- 
tentials that may redden the spectrum and reduce the 
string contribution to the CMB. 

We aim to formulate the model in such a way that the 
number of parameters remains small (e.g no more than 
7 parameters) and presume strings form. This can be 
achieved by the restriction to g < 0.1 and adding the 
potential term in the form of g 2 c 2 H S 2 . This way, the 
predictions remain independent of the gauge coupling g. 
We note that in presence of a negative mass-square term, 
there may also be models with g > 0.1, that are no longer 
disfavoured by the data. Unless one applies a restriction 
to more specific models, a systematic study of this region 
of parameter space would require additional free param- 
eters in the analysis and is therefore not performed here. 

In Fig. [2] we show 2D likelihoods for the non-minimal 
models using string spectrum A. In the F-term case, neg- 
ative values of c H lead to an increased red-tilt in ns, al- 



leviating the tension with data present for 



H 



0. The 



fit is slightly better than the 6 parameter model with 
A% 2 = 1.2. We find marginalized values of log 10 n = 
-2.28 ± 0.46, M = (0.346 ± 0.065) x 10 16 GeV and 



H 



= -0.027 ±0.011. 



For the D-term model the situation is more compli- 
cated, since for fixed Pr, and n the potential can admit 
two solutions, each with different M and parameters ns 
and Gfi. In the analysis wc treat these independently, 
running separate MCMC chains for each solution branch. 
The results of the analysis are shown in Fig. [5J The 
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FIG. 2: 2D likelihoods for non-minimal F- (top set of two) 
and D-term (bottom set of two) for string model A. In the 
D-term case the two sets of contours show the two solutions 
admitted by the potential, as explained in the text. 

solution for small k is essentially the same as the min- 
imal case. Here, M is higher and a non-zero c 2 H = 
does not effect the inflation dynamics. We find log 10 k = 
-4.45 ± 0.22 and M = (0.203 ± 0.032) x 10 16 GeV, with 
Ax 2 = -8.6. 

The other solution has a lower value of M, and a non- 
zero c 2 H flattens the potential near horizon exit. One gets 
sufficient c-foldings when the inflaton is placed close to 
the top of the potential and the resulting ng is redder, 
with the same slightly improved A% 2 = 1.2 as in the 
F-term model. For this branch larger values of k are al- 



lowed, since the lower M leads to a smaller string tension. 
We find log 10 k = -3.78 ± 0.44, M = (0.116 ± 0.071) x 
10 16 GeV and - log 10 (l0 18 c^/GcV 2 ) = 1.63 ± 1.50. We 
note that the upper bound on k is related to the upper 
bound onG/*. 

Conclusions: Minimal models of hybrid inflation rely 
only on two parameters (k and M), and they do not suf- 
fer from the potential problem of trans-Planckian values 
of the inflaton field. Therefore they are predictive and it 
is possible to derive parametric constraints from cosmo- 
logical data. In Refs. ^ has been found that a close 

to scale invariant spectrum of adiabatic perturbations is 
in agreement with CMB data, when a sizable string com- 
ponent is present. In particular, minimal models of hy- 
brid inflation which predict these features appeared to 
be viable in that light, a finding that may open inter- 
esting directions in model building [l6| . Using recent 
CMB data, SDSS and measurements of the deuterium 
abundance leads us to a revison of this picture: Mini- 
mal D-term models are ruled out at a level that can be 
estimated to be « 4a (the same is true for the HZ spec- 
trum) and the tension between minimal F-term models 
and the data is increasing. Non-minimal scenarios in- 
clude the negative mass-square term models discussed 
here and also the possibility of adding a negative-sign 
tadpole term [2l|. This way, one can entirely remove 
tension with the data, but at the same time, also the 
very distinct predictions of minimal hybrid inflation for 
the range of the spectral index are lost. 
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